A diameter— bandwidth product limitation of isolated-object cloaking 
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We show that cloaking of isolated objects is subject to a diameter-bandwidth product limitation: 
as the size of the object increases, the bandwidth of good (small cross-section) cloaking decreases 
inversely with the diameter, as a consequence of causality constraints even for perfect fabrication 
and materials with negligible absorption. This generalizes a previous result that perfect cloaking of 
isolated objects over a nonzero bandwidth violates causality. Furthermore, we demonstrate broader 
causality-based scaling limitations on any bandwidth-averaged cloaking cross-section, using complex 
analysis and the optical theorem to transform the frequency-averaged problem into a single scattering 
problem with transformed materials. 

PACS numbers: 



I. INTRODUCTION 

In this work, we extend the result that perfect isolated- 
object cloaking is impossible over a nonzero band- 
width pj, [^ to show that even imperfect cloaking of iso- 
lated objects necessarily has a bandwidth that decreases 
with the size of the object. More generally, we show that 
the cloaking efficiency (the ratio of total cross section 
to geometric cross section) must worsen proportional to 
the diameter when averaged over any finite bandwidth 
for cloaking of isolated objects in air (or any medium of 
negligible loss and dispersion). Unlike our previous proof 
that sensitivity to imperfections worsens with diameter 
at individual frequencies [^ , the result in this paper holds 
even for perfect fabrication and for materials with negli- 
gible absorption over the desired bandwidth. Our proof 
involves some unusual mathematical techniques. First, 
we equate the frequency-averaged problem to a single 
scattering problem at a complex frequency with the help 
of a version of the optical theorem. Second, we map 
the complex-frequency problem to an equivalent real- 
frequency problem with transformed materials. We can 
then use the fact that Ime and Im fi are > for any 
physical causal material at Re to and Im w > [J] to ana- 
lyze an "effective" absorption loss in a manner similar to 
our previous work |3|. Finally, in Sec. lIIII we numerically 
verify this scaling in an example of a spherical cloak that 
is a perfect Pendry cloak at one frequency but has causal 
dispersion. 

The idea of transformation-based invisibility cloaks 
was proposed in 2006 [l|, a fascinating idea that was 
followed by many theoretical works [5l-[28| and several ex- 
perimental demonstrations of cloaking of isolated objects 
at one frequency J29li32| . However, as pointed out by 
Pendry using a speed-of-light argument [1| and by Miller 
with a more formal approach [^, perfect cloaking of an 
isolated object in vacuum over a non-zero bandwidth is 
impossible due to causality. This severe limitation helped 
inspire the idea of ground-plane cloaking p33l - [37| and its 



experimental demonstrations J38M47| , which circumvents 
such causality limitations (although it is still subject to 
other practical scaling difficulties [3,|48[). However, the 
Pendry and Miller proofs say little about imperfect cloak- 
ing (a nonzero but small scattering cross-section). By 
continuity, if near-perfect cloaking is attained at a single 
frequency, "good" cloaking (total cross-section bounded 
by some given fraction of the geometric cross-section) 
must persist over some finite bandwidth. We show that 
causality imposes an even stronger constraint than for- 
bidding perfect cloaking over a finite bandwidth: for a 
bandwidth-limited cloak, we show that causality con- 
straints imply that the bandwidth of good cloaking scales 
inversely with the object diameter. (The impossibility of 
perfect cloaking over a finite bandwidth follows from our 
results as a special case.) Moreover, our proof holds for 
isolated objects in any transparent medium (negligible 
loss and dispersion), not just in vacuum. 

The key to our proof in Sec. |TT] is the transformation 
of a frequency-averaged scattering problem (weighted by 
a Lorentzian window for convenience), which is hard to 
analyze, to a single complex-frequency scattering prob- 
lem that is easy to analyze by relating the complex fre- 
quency to equivalent complex materials. In order to make 
this transformation, we rely on the optical theorem |49| . 
which relates the total cross section to the imaginary part 
of a forward scattering amplitude, since the latter is a 
causal linear response and hence an analytic function in 
the upper- half complex- frequency plane [J]. A review 
of the optical theorem and its history and applications 
can be found in [50]. The use of the optical theorem, 
and similar relationships based on conservation of en- 
ergy, to apply complex analysis to scattering problems 



is most common in quantum field theory 5l|. Similar 



in spirit to this paper, the optical theorem is also used 
in the "ITEP sum rules" of quantum chromodynamics 
to relate integrals of scattering cross sections multiplied 
by Lorentzian windows (or powers thereof) to scattering 
amplitudes at single points in momentum space via con- 



tour integration (5l|. On another related note, contour 
integration of the imaginary part of Green's functions 
(and other scattering ampUtudes) is also a key technique 
for computing Casimir interactions in quantum field the- 
ory |52| . where the relationship between the imaginary 
part of the Green's function and the field fluctuation 
statistics (the fluctuation-dissipation theorem) is, like 
the optical theorem, derived from energy-conservation 
considerations ISSll. 
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II. DERIVATION 

A transformation-based cloak, as depicted in Fig. f, 
works by mapping an object (in "physical" space X) to 
a single point (in "virtual space" X') and the cloak re- 
gion (volume Vc) to empty space (volume Vj,'), via a co- 
ordinate transformation that is the identity outside the 
cloak. If the ambient materials are Sa and /ia, then 
an ideal cloak is obtained by constructing the materi- 
als e — JSaJ' / dct J and /^ = J\iaJ' l^&'^J in the 
cloak, where J is the Jacobian matrix of the transfor- 
mation. In our previous work [3| we showed that the 
cloaking problem at one frequency becomes increasingly 
difhcult as the size of the object being cloaked increases. 
In particular, we proved that the imperfections due to 
absorption losses and random fabrication disorder must 
decrease asymptotically with the object diameter in or- 
der to maintain "good" cloaking performance: for the 
total (scattering -I- absorption) cross-section to be less 
than a given fraction / of a geometric cross section Sg. 
To prove these results, we assumed bounds on the attain- 
able refractive index contrast in the cloak, h < n/ua < B 
{n = y/ejl), which we showed to be equivalent to bounds 
on the singular values of J^. Using these bounds, we were 
then able to bound the field amplitudes in the cloak in 
terms of the field amplitudes in the virtual space (which 
are a constant for incident planewaves); such bounds, in 
turn, impose bounds on the allowed absorption losses and 
random disorder. In particular, we showed in the case of 
absorption losses that A Im e scales proportionally to yr , 
a measure of the inverse diameter. 

We now wish to analyze the effect of other imperfec- 
tions, especially inevitable material dispersion, on the 
cloaking performance, even in the idealized case of per- 
fect fabrication and negligible absorption. In particular, 
we suppose that one is interested in the average total 
cross section ctot over a bandwidth Aw around an oper- 
ating frequency Wop, which (for reasons described below) 
is convenient to define via a Lorentzian averaging weight 
as: 



FIG. 1: Schematic of transformation-based isolated-object 
cloak, which works by mapping the object and cloak in phys- 
ical space X to a virtual space X' in which the object is 
a single point. The transformation laws of Maxwell's equa- 
tions [l|? ] then dictate the required e and /i materials in the 
cloak volume to produce equivalent solutions in X and X' . 



and simply rescaling the whole system, we show that the 
cloaking problem becomes increasingly difficult, in the 
sense that {<Jtot)Au> I Sg ~ d (where ~ means proportional 
to), with Sg being a geometric cross sectional area of 
the object. (This scaling breaks down when CTtot/sg is 
no longer small, i.e. when it is no longer an approxi- 
mate cloak.) The only assumptions are that the ambient 
medium £a and /Iq is approximately lossless and disper- 
sionlcss over the bandwidth of interest (e.g. for cloaking 
in air) and that the attainable refractive index contrast 
(eigenvalues of ye^i/eaPa) is < B for some finite bound 
B (similar to our previous work [3|). 

Our analysis in the following subsections constitutes 
the following steps. First, in Sec. Ill A[ using complex 
analysis combined with the optical theorem, we relate 
the frequency average of Eq. ([T]) to a single scattering 
problem at a complex frequency so that we no longer 
need to consider the cross section at many frequencies at 
once. Second, in order to understand the precise meaning 
of this complex- frequency scattering problem, in Sec. Ill Bl 
we derive a variant of the optical theorem that is partic- 
ularly easy to analyze. Third, in Sec. IllCl we relate this 
complex-frequency scattering problem to an equivalent 
scattering problem at a real frequency with transformed 
complex materials. Fourth, in Sec. IIIDI we use analysis 
similar to our proof in [3| to show the "losses" introduced 
by these effective complex materials must scale with di- 
ameter, and hence mean crtot must scale with diameter. 
Finally, for the special case of a bandwidth-limited cloak 
(a cloak that is very good at one frequency but spoiled at 
other frequencies by material dispersion), we show that 
the bandwidth must narrow inversely with diameter in 
Secimi 
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As we increase the diameter d of the object and cloak, 
keeping the materials and transformation mapping fixed 



A. Frequency average of the scattering 
cross-section 

By the optical theorem [40|, (Jtot(w) = Im/(w) where 
/ is a forward-scattering amplitude (at least for the case 



of an incident planewave, but a generalization is given in 
the next section for any incident field), and so 



{crtot)ALj = Im / f{u) 



Alj/tt 



{uj — Wop)^ + Aa;2 



duj. (2) 



Because / is a causal linear response, it is analytic for 
Ima; > 0, and therefore we can perform a contour inte- 
gration to obtain (crtot)Acij in terms of the residue of the 
single pole of the Lorentzian in the upper-half complex 
plane: 



(o-tot>Aaj = Im/(wop + zAw). 



(3) 



This corresponds to a single scattering problem at a com- 
plex frequency, analyzed in more detail below. 



B. The optical theorem and analytic continuation 
to complex uj 

Consider any finite- volume scatterer (in a linear time- 
invariant system), described by some change Ae and A/i 
in the permittivity and the permeability compared to the 
ambient medium. Let the incident field (in the absence 
of the scatterer) be described by a six-component vector 



field '011, 



Hi 



In the presence of the scatterer. 



this is modified to a new total field -0 = ipinc + V'scat , the 
sum of the incident and scattered fields. The charge per- 
turbations in the scatterer are described by bound elec- 
tric and magnetic polarization currents AJ = — zwAeE 
and AK = — zwA/iH, respectively, which can be coni- 

AJ 



bincd into a six-component current field ^ 



AK 



Abstractly, we can write ^ = ^"^inc for some linear op- 
erator A relating the incident fields to the induced cur- 
rents, and causality (currents come after fields) implies 
that A is an analytic function in the upper-half complex- 
uj plane [J]. 

Physically, the scattered field is the field produced by 
these oscillating induced currents ^ in the scatterer, and 
the interactions of the currents and fields provide a sim- 
ple way to characterize the absorbed and scattered pow- 
ers. The total absorbed power Pabs, assuming an ambi- 
ent medium with negligible dissipation, equals the total 
time-average {tp) incoming Poynting flux: 



Pabs = ^ Rc y"(E* • AJ + H* • AK) 

= ^ Re (0,0 

= 2 Re (0inc, + ^ RC (0scat, : 



(4) 



where we have defined the inner product (• • • ,•••). Sim- 
ilarly, the scattered power is the work done by the cur- 
rents on the scattered field: — ^ Re (V'scatiO- Therefore, 



the total (absorbed + scattered) power is: 

-'tot -* abs T Jgcat 



= Q Re (0inc, + ^ Re (V'scat, j - ^ Re (V'scat, Cl 

= 2Re(-0inc,O 

= -Rc(V'inc, A0inc) 



Im/(w), 



(5) 



where we have defined /(w) = | (^inc, ^V^inc)- The in- 
terpretation of / as a forward-scattering amplitude when 
"01110 is a planewave is actually irrelevant to our proof, but 
it follows from the fact that / in that case is simply the 
forward-planewave Fourier component of the field r^ iS,/uj 
corresponding to the currents ^. 

The key question, for our purposes, is to understand 
what / looks like at a complex frequency. Suppose we 
have homogeneous, lossless ambient medium with an in- 
cident planewave ip'mc = '0oe*'^^' "^ propagating in the -l-a; 
direction for a constant amplitude ^o- Then 



/(^) 



-00,6' 



-iujx I c 



A{ijj)e 



-^-iujx/ c 



-00 



(6) 



moving all the x dependence to the right-hand side of the 
inner product. At a complex frequency w — )■ a; + i7 in the 
upper-half plane (7 > 0), A{uj + i"f) is analytic by causal- 
ity, while the planewave terms are analytic everywhere 
(and become exponentially decaying or growing waves in 
space at complex a;). So, / at a complex frequency rep- 
resents an overlap with an exponentially growing field of 
the currents produced (via A) by an exponentially de- 
caying source. In the next section, we clarify this picture 
further by relating the scattering operator A at a complex 
frequency to scattering at a real frequency with complex 
materials. 



In order to perform contour integrations of f{iS) mul- 
tiplied by a Lorentzian, it is not enough for / to be an- 
alytic, however: it must also be bounded as |w| — > 00 
so that we can close the contour above. This fact is 
used extensively in quantum field theory [5l|. Intu- 
itively, this occurs because the the exponential growth of 
g-l-7x/c cancels the exponential decay of e"''^/'^. Math- 
ematically, as |a;| — >■ 00 the Ae and A/i must vanish 
(since all susceptibilities vanish in the limit of infinite fre- 
quency [1, |40|), and so A simplifies: to lowest order for 
weak scatterers (i.e., in the first Born approximation), 
— zwAeEjnc 

-ZwA/lHinc 

ponential factors exactly cancel. We will use a similar 
procedure, below, to analyze the effects of small imper- 
fections in the cloak due to material dispersion. 



C ~ ^0i, 



in which case the ex- 



C. From complex frequencies to complex materials 



1. Irrelevance of artificial absorption 



Combining the previous two seetions, we can now re- 
late the frequency-averaged scattering cross-section, for 
an incident planewave (say in the x direction to the so- 
lution of a single scattering problem at a single complex 
frequency: 



(o-tot>A(^ = 7:Rc('i/'o,e 



^ — iujx I c 



A{uj). 



^-\-iujx I c 



A 



(7) 



evaluated at a; = Wop + iAw. The A(a;)e+*'^'^/'^Vo repre- 
sents the induced currents ^, for the materials evaluated 
at the complex oj, in response to an exponentially decay- 
ing incident planewave (multiplied by a constant ampli- 
tude -00 )■ The interpretation of this problem is simplified 
by the fact that a complex frequency is mathematically 
equivalent to using modified materials at a real frequency 
Wop- In particular, consider Maxwell's equations in the 
frequency domain: 



V X E = -iuj^iU - K, 

V X H = +iujeE + J, 



(8) 
(9) 



where J and K are (free) electric and magnetic cur- 
rent densities, respectively. For our complex w, the 
iweE term becomes i{ujop + iAuj)e{ijJop + JAu;)E ~ 



e{oJop + iAio).{l + i^) 



WopE, 



which looks exactly 

like the term for a real frequency ujop at a modified per- 
mittivity i 



The first change in the materials is the multiplica- 
tion of e and /i by 1 + i-^. But this change does 
not hurt the cloaking performance because it is done 
uniformly everywhere in space. More explicitly, if e = 
J^aJ"^ I det jj is a valid transformation-based cloak, so 
is £.(1 + ^) = j£ah + ^)J^ I det jr, except that this 
is now a transformation-based cloak for a lossy ambient 
medium £a.(l + i^), and similarly for ^. Therefore this 
change in the materials leaves the scattering cross-section 
invariant. 



2. Impact of material dispersion 



The second change in the materials is that we need to 
evaluate e and fi at the complex frequency Wop -I- iAuj in- 
stead of Wop, and here the presence of material dispersion 
(unavoidable for any material other than vacuum) acts 
to spoil the cloak. In particular, as reviewed in the Ap- 
pendix, causality and other fundamental principles imply 
that Ime and Im/z are both strictly positive at cUop + iAuj 
for Ao; > Q, corresponding to an unavoidable ad- 
ditional absorption in the complex-frequency scattering 
problem. Unlike the artificial absorption in the previous 
section, this is an absorption defect introduced only in 
the cloak: e (both Im and Re) differs from the cloaking 
transformation of the ambient medium by 



e(wop) = e(wop + iAuj). (l + i 



Auj 



(10) 



Similarly for fi ^ fl. 

Thus, operating at a complex frequency is equivalent to 
two changes in the materials. First, the multiplication by 
l+iAtu/oJop corresponds to an effective absorption added 
throughout all space (including the ambient medium). 
Second we evaluate e and fj, at Wop + iAuj rather than at 
the real frequency ojop, which will change their values in 
the presence of material dispersion. 



D. Consequences of dispersion on 
frequency-averaged scattering 

So far. we have shown that the problem of finding the 
frequency-averaged scattering cross section is equivalent 
to solving a certain scattering problem at a single com- 
plex frequency, which in turn is equivalent to a scattering 
problem at a single real frequency Wop with modified com- 
plex materials e and fl. In this section, we analyze the 
consequences of those effective material modifications, di- 
vided into two separate material changes as discussed 
above, for cloaking performance. 



Ae = [£(a;op + lAuj) - £(a;op)] ( 1 + i-^ 



and similarly for /i. The key assumptions here are that 
the ambient medium has no dispersion over the given 
bandwidth, so that the ideal cloaking transformation at 
the complex uj is given by the second term in A£, and that 
the ambient medium is lossless, so that e{uJop) is real and 
ImA£ > 0. More generally, it is sufficient for the disper- 
sion and loss of the ambient medium to be small enough, 
compared to the dispersion of the cloak materials, such 
that ImA£ is > 0. 

We can analyze the consequences of this effective ab- 
sorption imperfection similarly to our previous work |3| 
It is convenient to first transform the imperfections to 
virtual space (as shown in Fig. [ij where the scatter- 
ing problem is easier, because in the absence of Ae and 
A/i we have a planewave in a homogeneous medium. In 
particular, given an upper bound B on the index con- 
trast as discussed above, one obtains Ae' > Ae/B in 
virtual space, and similarly for /i [3|. Furthermore we 
are only interested in the regime in which we have a 
good cloak — once the imperfections become so large as to 
make the cloak useless, all of the scaling relations break 
down and the problem is no longer interesting. This is 
the regime in which Ae and A/x are small, and there- 
fore virtual space at complex ui is equivalent to a ho- 



notation of Sec. HTbI C = ^Vii 



mogcncous medium with small imperfections and per- 
turbative methods are applicable. In particular, the 
lowest-order scattering current (in virtual space) is sim- 
ply J' = —iujopAe'Einc, and similarly for K'. In the 

-iojA/i'Hin 

Substituting this into Eq. ([7]), we find that the exponen- 
tial factors e^'^^'^^/'^ exponential factors exactly cancel, 
leaving: 

(atot)Ac.«^/ [lEol'lmAe' + IHol'lmA/i']. (11) 

This has two main consequences. First, meancrtot > 
for Aw > since Im Ae and Im A/z are strictly positive as 
noted above: even if the cloak is a perfect cloak at Wop, it 
is imperfect when averaged over any non-zero bandwidth. 
This is, therefore, an alternative proof of the results of 
Pcndry [1| and Miller [2| that cloaking of isolated objects 
over a non-zero bandwidth is impossible for physical, 
causal materials. Second, exactly as we showed for other 
imperfections in previous work [3[, it immediately fol- 
lows that ((Ttot)Ai^ grows ^ V^ ^ Vc ^ Vo and hence the 
frequency-averaged cloaking efficiency {(Jtot) Alo/ Sg scales 
proportionally to a mean diameter Vo/sg. This is the 
central result of our proof, but to better understand its 
consequences we consider some special cases in the next 
section. 



E. Scaling of cloaking bandwidth \vith diameter 

As shown in the previous section, the fractional cross- 
section mcancTtot/sg, averaged over a bandwidth Aw 



around 



-"op J 



must scale proportional to the diameter d. 



at least as long as {(Jtot)Auj/sg ^ 1 (i.e. until cloaking 
breaks down completely). There are two possible sources 
of this linear scaling, depending on whether (fTtot)Ati; is 
limited by the cross-section at oJop (due to imperfections 
in the cloak at the design frequency) or by the bandwidth 
of a dip in the cross-section around Wop (due to mate- 
rial dispersion degrading a near-perfect single-frequency 
cloak). In the former case, the physical mechanism is 
simply that the losses due to imperfections at Wop scale 
with diameter, as we already proved in Ref. y. In the 
latter case, however, it leads to a new prediction: in a 
bandwidth-limited cloak, the bandwidth must narrow as 
the object diameter increases. In fact, we show in this 
section that the bandwidth gcnerically narrows inversely 
with the diameter in this case. 

First, let us consider the bandwidth scaling from 
generic dimensional considerations. Because of material 
dispersion, one expects good cloaking to only be possible 
in some limited bandwidth ~' F around some design fre- 
quency Wop, in which case a small meantrtot/sg requires 
Aw <C F. If we Taylor-expand {<Jtot) AujI Sg in Aw/F <C 1, 



we would generically expect an expansion of the form: 

{(Ttot) Au/sg = crtot('^op)/sg + CAw/F + 0[(Aw/F)^] 

for some coefficient C, where C is generically > if Wop 
is chosen to be a minimum of crtot • For sufficiently small 
bandwidths Aw ^ rcrtot(wop)/sg, this is dominated by 
the scattering at Wop, which scales linearly with diameter 
in the presence of imperfections as shown in Ref. y. On 
the other hand, for a good single- frequency cloak at Wop, 
there is a regime F(Ttot('^op)/sg ^ Aw <C F where the 
second term dominates, i.e. where material dispersion is 
the limiting factor. The CAw/F term must therefore also 
scale linearly with the diameter in order for (crtot) Aw /sg 
to scale linearly for such Aw, and hence we can conclude 
that the bandwidth F of the cloak generally scales in- 
versely with diameter. 

This is best illustrated by a simple example. Consider 
the case where atot/sg achieves a minimum / ^ 1 at Wop 
with an approximately Lorentzian lineshape of width F, 
going to 1 at frequencies far from Wop: 



O-tot(w) 



(W - Wop)2 + F2 



(I-/)- (12) 



The integral of Eq. ([T]) for this crtot(w) can be evaluated 
analytically to obtain 



(o'tot)Aw = 1 — 



1-/ 



1 + Aw/F 



/+(l-/)Aw/F+0[(Aw/F)2 



Since this must scale linearly with diameter for any / <C 1 
and any Aw <C F, it follows that both / and 1/F scale 
linearly with diameter until cloaking breaks down. 



III. NUMERICAL EXAMPLE 

To illustrate and validate our predictions of this scaling 
of cloaking bandwidth, we performed explicit numerical 
calculations of the scaling for an example bandwidth- 
limited spherical-cloaking problem (near-perfect at one 
frequency, but degraded by causal dispersion at other 
frequencies) in vacuum (with nondimensionalized units 
Ea^eo^l and ^„ = /^q = !)■ 

At an operating frequency Wop, we use an exact 
"Pendry" cloak [Ij : a sphere of radius Ri is surrounded 
by a cloak of radius R2 > Ri that is linearly mapped to 
an empty sphere {R'l — 0), resulting in materials: 



Sri^op) — Mr('^op) 



R2 — Ri 
R2 



Ri 



£e(wop) = £<p{i^op) = fJ-e i^op) = /^^(wop) 



R2 



i?2 — -Rl 



We fixed i?2 = 1.5i?i, so that the cloaking material pa- 
rameters are the same for all Ri , merely rescaled in space 




FIG. 2; Relative cross-section versus frequency for a spher- 
ical cloak designed to be a perfect Pendry cloak at Wop and 
showing the effects of material dispersion at other frequen- 
cies, computed by a spectral scattering-matrix method. As 
predicted, the cloaking bandwidth decreases linearly with the 
object radius, for three object radii relative to Aop ~ 27rc/ciJop. 



bandwidth around Wop , and this bandwidth becomes nar- 
rower as the diameter increases. Quantitatively, if we 
look at the bandwidth at a fixed cr/sg of about 1/4 its 
maximum, we find that the bandwidths for Ri = 2Aop 
and i?i = 4Aop are « 1/2.1 and 1/4.1 times the band- 
width for i?i = Aop, respectively, almost exactly the pre- 
dicted linear scaling. 

[Far from Wop, one expects ctot to approach twice the 
total geometric cross section, in the limit of a large scat- 
terer, because in this Hmit all of the incident light is 
scattered. The factor of two comes from the definition of 
scattering cross section [40| , in which the scattered power 
appears twice: once as the scattered waves propagating 
in other directions, and once as the "shadow" canceling 
the forward-propagating wave. Here, the total geomet- 
ric cross-section is that of the cloak, 7ri?|, so we expect 
cJtot/sg « 101og[2(7ri?|)/(7ri?f)] « 6.5 (dB) away from 
Wop for i?i — > oo, and this is roughly what is seen in 
Fig.H] 



as the object becomes larger or smaller. 

By construction, at Wop we have a perfect cloak, and 
the only limiting factor is the bandwidth: we use ideal- 
ized lossless materials but with causal dispersion relations 
that satisfy the Kramers-Kronig constraints. In particu- 
lar, we use a combination of two limiting cases: a plasma 
model (a limit of a Drude model as losses go to zero) and 
a limit of lossless Lorentzian resonance (corresponding 
to a polarization field described by a lossless harmonic 
oscillator) at a frequency uq. Combined with the pre- 
scribed Pendry values at Wop from above, this results in 
dispersion relations: 



^(r,w) = ^-Lr{r,uj) = 1- 



£e — £0 — Me — M0 



1- 



i?2 



i?2 ^ Rl 



Ri 



r- Rl 



w 



op 
2 ' 



i?2 -Rl l^i- w2 



The Lorentzian resonance frequency ujq can be chosen 
arbitrarily; we used wq = 2a;op. 

This geometry was then simulated using a spec- 
tral (spherical-harmonic expansion) scattering-matrix 
method as described in Ref. |54j. The continuously vary- 
ing anisotropic material parameters are approximated 
by a large number of piecewise-homogeneous isotropic 
layers [5J|. The total scattering cross-section was com- 
puted over a range of frequencies for i?i = Aop, 2Aop, 4Aop 
(where Aop = 27rc/wop); and is plotted in Fig.[2j The re- 
sults in Fig. [2] are converged with resolution (the number 
of spherical layers) to within a few percent accuracy. At 
Wop, the cloak should theoretically be perfect, but we ob- 
tain a small nonzero (J / Sg (< 10^"^) due to the discretiza- 
tion errors, which vanishes with increasing resolution. 

As expected, the material dispersion prevents this from 
being a good cloak except at frequencies in a narrow 



IV. CONCLUDING REMARKS 

{'^totlAuj In this work, we extended on our previous 
paper [3| to study the bandwidth limitation of isolated- 
object cloaking which is a key limiting factor for this type 
of cloaking [l|, ^ . Although it was known that perfect 
cloaking was impossible over a nonzero bandwidth, this 
result did not seem to exclude the possibility of imper- 
fect cloaking over a finite bandwidth. Indeed, imperfect 
finite-bandwidth cloaking is possible, but we have now 
shown that it is subject to a severe practical constraint: 
the bandwidth inevitably narrows proportional to the ob- 
ject diameter, given fixed materials. Although we cannot 
infer any hard upper bounds on the size or bandwidth of 
such cloaking without further information about the at- 
tainable materials, this result indicates a fundamental 
challenge in scaling up small experimental demonstra- 
tions to larger cloaks. 

The use of gain has been prop osed to compensate 
for loss problems in cloaking |28l ? ]. Although gain 
is necessarily nonlinear and can be detected by a suf- 
ficiently strong incident field, in the idealization of lin- 
ear gain then many of the techniques in this paper are 
complicated by the fact that a linear-gain resonance (the 
complex-conjugate of an absorption resonance) would be 
non-analytic (and have negative imaginary parts) in the 
upper-half complex-frequency plane. (The average crtot 
must also be replace by the average |crtot| or similar, since 
gain can produce a dtot < 0.) However, the bandwidth 
limitations described in this paper arise even for idealized 
materials with negligible dissipation loss, due to disper- 
sion in the real part of the e and ^ alone, in which case 
even idealized linear gain is inapplicable. (Furthermore, 
as pointed out in our previous work [3|], gain compensa- 
tion of absorption must become increasingly perfect as 
the object diameter increases, nor does it compensate for 
scattering from fabrication disorder.) 



Ground-plane cloaking does not suffer from any intrin- 
sic limitation on its bandwidth from causality, nor does 
our proof apply in that case. The reason our proof docs 
not apply (at least, in the present form) to ground-plane 
cloaking is the figure of merit is no longer crtot : a ground- 
plane is actually designed to reflect waves, albeit to re- 
flect them in a way that mimics the ground plane. On the 
other hand, our previous work j3|, i48j] showed that even 
ground-plane cloaks are increasingly sensitive to imper- 
fections as the size of the object increases. This suggests 
that ground-plane cloaks should also be increasingly sen- 
sitive to material dispersion, for cloaking over a finite 
bandwidth, as the object size increases, and a challenge 
for future work is to quantify (or disprove) this relation- 
ship. 

An alternative direction for cloaking theory (and ex- 
periment) is to consider relaxations of the cloaking prob- 
lem that might prove more practical. In particular, it 
would be valuable to make precise the intuition that the 
cloaking problem becomes easier if the incident waves are 
restricted (e.g. to plane waves from a certain range of an- 
gles) and/or the observer is limited (e.g. only scattered 
waves at certain angles are visible, or only amplitude but 
not phase can be detected, or sufficiently small time de- 
lays are undetectable), since this is arguably the situation 
in most experiments. For example, current "stealth" air- 
craft are designed in the radar regime mainly to reduce 
back-scattering only [55|- So, it is clear that a sufficiently 
relaxed cloaking problem is practical even for large ob- 
jects, and one interesting goal is to find the "weakest" 
relaxation that remains practical at useful scales. 
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Appendix 

Here, we review a known consequence of causality that 
is the key to our analysis above: for a passive medium, 
Ime(a;) > in the upper half plane Imw > as proved 
inRef.lj. A condensed proof of this fact is as follows: luie 
is analytic in w in the upper-half plane by causality, and 
is therefore a harmonic function in the upper-half plane, 
and so can obtains its minimum only on the boundary of 
its domain, except in the trivial case of vacuum where it 
is a constant function (Ime = 0). In particular, consider 
the upper-right quadrant of the complex-w plane. Along 
the positive real axis, Ime > for a passive material (in 
the absence of gain), even for idealized lossless materi- 
als. Along the positive imaginary axis, Ime = since 
e(— a;) = e(i^)* for real- valued physical fields [43 . As 
|aj| — > oo one must have Imw — >■ 0. Hence the minimum 
of Im e along the boundary of the upper-right quadrant 
is zero and it is strictly positive in the interior. 

For physical materials, Imw > along the positive 
real-w axis except at w = in order to satisfy the sec- 
ond law of thermodynamics, and this is the usual case 
in which the above statement is proved |J|. However, it 
is also interesting to consider the idealized limit of loss- 
less materials (such as a plasma model or a lossless res- 
onance), in order to study bandwidth-averaged cloaking 
with idealized materials. Our proof, above, works even in 
this case as long as one is a little cautious about the case 
of poles in e lying exactly on the real-w axis, in order to 
exclude the case where Ime diverges to — oo as the real 
axis is approached from above. In particular, we must 
restrict ourselves to idealized materials that are the limit 
of physical lossy materials as the losses go zero, so that 
we are taking the limit as poles approach the real axis 
from below. In this case, Ime in the upper-right quad- 
rants is the limit of a strictly positive quantity and hence 
cannot go to — cxd. 
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